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New soliton solutions for thick branes in 4 + 1 dimensions are considered in this article. 
In particular, brane models based on the sine-Gordon (SG), and scalar fields are 
investigated; in some cases Z 2 symmetry is broken. Besides, these soliton solutions are 
responsible for supporting and stabilizing the thick branes. In these models, the origin 
of the symmetry breaking resides in the fact that the modified scalar field potential may 
have non-degenerate vacuua and these non-degenerate vacuua determine the cosmological 
constant on both sides of the brane. At last, in order to explore the particle motion in the 
neighborhood of the brane, the geodesic equations along the fifth dimension are studied. 
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1. Introduction 

Since there is no known fundamental principle requiring spacetime to be (3 + 
1)— dimensional— , our observable universe might be a (3 + 1)—dimensional brane 
in a higher dimensional bulk— Most models suggest that there are one or more 
flat 3-branes embedded discontinuously in the ambient geometry—. Furthermore, 
almost all of the extra-dimensional models require the existence of scalar fields, for 
instance, to generate a domain-wall which localizes matter fields^. Scalar fields 
serve to stabilize the size of the compact extra dimensions^, and can support the 
modification of the Randall-Sundrum warped-space^ to a smoothed-out version— 
or to cut off the extra dimension at a singularity—!^. In addition to this, the non¬ 
linearity in the scalar field and, in particular, the existence of disconnected vacuua 
in the self-interaction of the scalar field lead to the appearance of a stable localized 
solution, which is a good motivation for building thick brane models^. It is fre¬ 
quently invoked to replace an infinitely thin brane with a thick one, by supposing 
a scalar field with soliton behavior. Based on this point of view, in this paper we 
investigate diverse models, namely, the sine-Gordon (SG), ip'^ and brane models 
which have broken Z 2 symmetry in some cases. Note that the Z 2 symmetry may be 
restored by a proper choice of model parameters. The origin of symmetry breaking 
in our models reside in the fact that the modified scalar field potential may have 
non-degenerate vacuua.^. The cosmological constant on both sides of the brane are 
determine by these vacuua. 
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We consider a thick brane, embedded in a five-dimensional (5D) bulk spacetime, 
modelled by the following action: 


S' = y d^xy/\g\ 




( 1 ) 


The simplest line element of the brane, embedded in the five-dimensional bulk 
space-time can be written as^: 


ds\ = gAsdx^dx^ 

= dw"^ + e^^{dx‘^ + dy^ + dz^ — dt^), 


( 2 ) 


where A is the warp factor which depends only on the 5D coordinate w. 

The 5D gravitational and scalar field equations take respectively the following 
forms 


2,A" -h = -nie-^^Too = - 


+ V{lp) 


= kIt^a = 




ff I A Af f 

Lp + AiA (p = 


dV{^) ^ 

dip 


( 3 ) 

( 4 ) 

( 5 ) 


where the prime denotes derivative with respect to w. 

In order to obtain the first-order equation, one can introduce the auxiliary func¬ 
tion W in such a way that: 


A' = -lw(ip), 

, _ 1 dW(ip) 

2 dp, ^ 

while V(ip) achieves the special form— 


V(ip) 


_ 1 / dW(ip) Y _ 1 


w(p,y. 


( 6 ) 


( 7 ) 


8 \ dp J 

Moreover, it may be instructive to calculate the geodesic equation along the fifth 
dimension in a thick brane, in order to investigate the particle motion near the 
brane^. To this end, we start with the geodesic equation: 




0 dx^ dx^ n 

+ ^ AB — -^— — 'J 


dr dr 
4 dx^ dx^ 
dr'^ dr dr 




d?s 


= 0 


w -k = 0 


which leads to 

w + clf{w) = 0, 

where ci is a constant of integration and the function f(w) is defined as 

f{w) = A'{w)e-^^^^\ 


( 8 ) 

(9) 


(10) 
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2 . Thick brane models 

First we consider the Sine-Gordon (SG) system, where the self-interaction potential 
of this model reads: 


cos{b<f)] 


( 11 ) 


where a and b are free parameters of the model. The SG system has the following 
exact static kink solution—: 



( 12 ) 


which is plotted in FiglTJa), for various values of parameters a and b, which corre¬ 
spond to branes with different thickness. When considered as the brane potential, 
however, this potential should be modified to become consistent with the Einstein 
equations. The corresponding potential for this model is given by: 



(13) 


which is depicted in Figl^Ka). Notice that this potential has two series of non¬ 
degenerate vacuua, as in the DSG (double sine-Gordon) system potential—. How¬ 
ever, in the limit oib ^ a these vacuua tend to the same value (become degenerate), 
such as the potentials used in—lii. 

In this case the energy density is localized at the brane and the thickness of the 
latter is A = It can be seen that there is no singularity in the Ricci scalar 

and/or Kretschmann scalar. Moreover, in the limits of ic —> ±oo, the Ricci scalar 
and all the components of the Einstein in the right and left sides of the brane are 
different (Figl^Ia)). However, in the limit of w —>■ 0, the Einstein tensor is given 
by G/ = I (for ^ = v = 0,1,2, 3), and one can interpret it as the 

cosmological constant on the brane, i.e., G® oc A5®. 

Therefore, we have a broken Z 2 -symmetry in the bulk, as the two sides of the 
brane differ completely. On the right and in the limit of ic —>■ +oo, the Einstein 
tensor and consequently the cosmological constant of the bulk vanish, so the bulk is 
asymptotically Minkowski. However, on the other side of the brane, these quantities 
are nonzero and equal to the constant value 512a/(36^), and as a result the bulk 
would be de Sitter. Besides, by checking the geodesic equation of a test particle 
moving only in the direction of the extra dimension, the confining effect of the scalar 
field is proved. 

The second brane model is the (/s'^-based model, for which we have>^: 



(14) 


where a and /3 are constants. The kink solution reads 


^p{w) = atanh{Pw), 


( 15 ) 













V{<f) 
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Fig. 1. Soliton solutions as a function of the fifth dimension for the following models: (a) SG for 
6=1, (b) (j)'^ for Q = 1 and (c) for a = 1 systems. Dashed, dotted-dashed, and continuous 
curves correspond to solitons with decreasing brane thickness. In the limit of an infinite a//3 
parameter, the soliton approaches the step function. 


which is depicted in FiglTJb). The potential is obtained as 

2 


yi^) = 


1 - 


‘p 


^ 2 2 /d2 


3- 


(16) 


which is plotted in Figl2](b) and the brane thickness becomes A = As for 

the previous SG model, we determine the limits w —>■ ±oo for all the components 
of the Einstein tensor components, which are the same at different sides and given 
by (FigEKb)) and in the limit of w —>■ 0 the Einstein tensor components 

(for ^ = 1 / = 0,1,2, 3) takes the form Therefore, in this model the 

cosmological constant on the brane A would be For this model a test 

particle oscillate around the brane by 11 = -y/1 




P P P 


Fig. 2. The plots depict the modified soliton potential as a function of the fifth dimension for (a) 
SG with a = b = 10, (b) with a = 0 = 1 and (c) with a = 0 = 1 systems. The potential of 
the SG and systems have non-degenerate vacuua. In contrast, the (0^ potential has degenerate 
vacuua and this leads to stable, topological solitonic brane. 
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The last model is based on the soliton. For this model, we have the following 
potential: 

«2 




(17) 


where a and /3 are constant(as in the model) and as a result, the kink solution 
is given by—: 

a 


(i){w) = 


(18) 


The kink solution is depicted in FiglUc). Thus, for the system the self-interaction 
potential takes the form (FigEJc)): 

vm = 7 — 


(19) 

The thickness of this brane is given by A = . Moreover, for this model 

and in the limit of w —> ±oo the mixed Einstein tensor components reduce to 
and zero respectively see Figl^Kc). and in the limit of w —>■ 0, the Einstein 
tensor is ~ l) 



WWW 


Fig. 3. The plots depict the Einstein tensor component Gq for the (a) SG with a = b = 1, (b) (p^ 
with a = 0 = 1 and (c) (p^ with a = 0 = 1 systems. Note that this quantity approaches different 
constant values for the SG and systems, while the system is Z 2 -symmetric. 


3. Conclusion 

In this work, we obtained exact thick brane models inspired by well-known nonlinear 
systems, namely, the sine-Gordon (SG), and models. The confining effect of 
the scalar field in all these three models were confirmed by examining the geodesic 
equation for a test particle moving normal to the brane. In particular, it turns out 
that the modified potential for the SG system resembles that of the double sine- 
Gordon {DSG) system, while those of and became and respectively. 
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In the case of the SG model, the resulting brane does not have Z 2 symmetry, 
in general, where the center of the brane may be displaced from w = 0 and the 
potential will not be an odd function of w in general. However, by a suitable choice 
of the model parameters it is possible to make the vacuua of the effective potential 
degenerate, in which case the Z 2 symmetry is restored. In the case of the model, 
however, we could not restore this symmetry via re-parametrization. 
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